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Abstract 

A spatially-discrete sine-Gordon system with some novel features is described. There is a 
topological or Bogomol'nyi lower bound on the energy of a kink, and an explicit static kink 
which saturates this bound. There is no Peierls potential barrier, and consequently the motion 
of a kink is simpler, especially at low speeds. At higher speeds, it radiates and slows down. 

AMS Classification numbers 58F, 70F. 

1 Introduction 

There are many nonlinear systems, in one spatial dimension, which admit topologically-stable kink 
solutions: for example, the sine-Gordon and phi-four systems. They have many physical applications. 
For applications in, say, condensed-matter physics or biophysics, an accurate model should take the 
discreteness of space into account — in other words, the kinks live on a one-dimensional lattice rather 
than in a continuum. A discrete version of the sine-Gordon system which has been extensively studied 
is the Frenkel-Kontorova model, in which the partial derivative d/dx of the continuum system is 
replaced by the forward difference on the lattice. In this case, much of the topological character of 
the model is lost, and there is no known explicit kink solution. 

Let us briefly recall two of the features of kink motion in the Frenkel-Kontorova model [1~5]. 
First, the energy of a kink depends on its position in relation to the lattice: a static kink located 
halfway between two lattice sites has a lower energy than one located exactly on a lattice site. The 
difference between these two energies is the so-called Peierls-Nabarro barrier. So if a moving kink 
does not have enough kinetic energy, then it gets trapped between two adjacent lattice sites; it 
oscillates, and emits radiation, and cannot escape. Secondly, a kink which starts off with enough 
speed to surmount the Peierls-Nabarro barrier, and move along the lattice, will lose energy (as 
radiation) and slow down, until eventually it too becomes trapped. 

In this note, we wish to describe an alternative discrete sine-Gordon system, in which the Peierls- 
Nabarro barrier is eliminated. It maintains an important feature of the continuum model, namely a 
topological lower bound ( "Bogomol'nyi bound" ) on the energy of a kink. And it admits an explicit 
static kink solution on the lattice, which saturates this lower bound. A moving kink will still radiate 
and slow down, but it is never trapped, no matter how slowly it moves. So this lattice sine-Gordon 
model is quite different from the usual one described above. 

In the next section we shall describe the model. Subsequent sections study kink motion, for low 
speed and for (moderately) high speed. In each case, we make approximate analytic predictions of 
the kink behaviour, and compare these with numerical simulations. 
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2 A topological discrete sine-Gordon system 



Let us begin with a very brief review of the situation for the continuous sine-Gordon equation 
'^xx ~ 'Ptt — sin (p. It is convenient to use the dependent variable tp — ^tp instead of ip. The potential 
energy of the field is 

1 /"^ 

EP^7 {^l+Bm'^)dx (1) 



(this is normalized so that a single static kink has unit energy). The kink boundary condition is 
— *■ as a; — >• — oo, xp —> it as x —>■ oo. The standard Bogomol'nyi argument |^ is 

1 f°° 



1 f°° 

Ep + - da:{cosi;) dx 

^ J -oo 

Ep - 1, (2) 



where kink boundary conditions are imposed. So the energy Ep is bounded below by 1; and Ep 
equals 1 if and only if — smip, the solution of which is the static kink 

V'(a;) = 2arctanexp(x — a;o)- (3) 

From this point on, x becomes a discrete variable, with lattice spacing h. The subscript + 
denotes forward shift, ie. f+{x) = f{x + h); so, for example, the forward difference is given by 
A/ = h~^{f+ — /). To obtain a lattice version of the Bogomol'nyi bound, we may begin with the 
same function cos-^ as appears in (2), and "reconstruct" the inequality. The first step is to choose 
a factorization 

Acosip ^-DF, (4) 

where D —f ip^ and F —>■ sin?/; in the continuum limit /i — > 0. Then define the potential energy of 
the lattice sine-Gordon field to be 

, oo 

^^-i E (D'+F'). (5) 



a; — — oo 



It follows, just as in the continuum case, that (with kink boundary conditions) Ep is bounded below 
by 1; and the minimum is attained if and only if D = F. 

There is a choice involved in (4), and we are free to make whatever choice we like. The most 
natural one seems to be 

2 1 
D = -sin-(V'+ -V'), 

F = sini(^++V^). (6) 

Substituting these into (5) gives 

h ^ 4: \ 1 

which, of course, reduces to (1) in the continuum limit. 

Let us summarize what we have so far. The real-valued field ilj{t,x) depends on the continuous 
variable t and the discrete variable x. Its potential energy Ep is defined by the expression (7); for 
the kinetic energy, we may for example make the simple choice 

J oo 

a; — — OO 
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where ip = dtp / dt. The boundary condition on is that it should tend to an integer multiple of tt 
(independent of time) as a; — > — oo or x — > +cxd; this guarantees finite energy. A 1-kink configuration 
has the boundary values — > as x — > — oo and ip ^ it ds x ^ +00. For such fields, the total 
energy E = Ep + Ek is bounded below by 1; and this lower bound is attained if and only if = 0, 
and the two expressions in (6) are equal. 

This latter condition, namely D = F, is called the Bogomol'nyi equation ||]. It is a first-order 
difference equation, whose solutions (for kink boundary conditions) minimize the potential energy. 
Hence these (static) solutions are also solutions of the Euler-Lagrange equations 

2 dEp 

since dEp /dip = at a minimum. This is a general feature of the Bogomol'nyi argument: one gets 
first-order equations whose solutions are also static solutions of the second-order equations of motion. 
Furthermore (and this is the most important feature), solutions of the Bogomol'nyi equations have 
an energy which is at its topological minimum value. 

In the expression (7) for i?p, the first term in the summand represents an attractive force 
(nonlinear Hooke's law) which couples nearest neighbours, and the second term a substrate potential 
depending on the average of pairs of nearest neighbours. The strengths of these two forces have, 
in effect, been normalized by scaling, and h may be regarded as a dimensionless parameter in the 
model. 

The Bogomol'nyi equation D = F may also be written as 

V+ 2 + /1 V 
tan— — = -tan—, (10) 

from which one sees that the parameter h should be less than 2 if one is to obtain a well-behaved 
solution. In fact, the solution of (10) can be written down exphcitly: given that < ft. < 2, it is 

■i/' = 2 arctanexpa(x — 6), (11) 

where 

a = — arctanh— , (12) 

h 2 ^ ^ 

and where b is an arbitrary real constant. This is a static lattice kink solution. Because it corresponds 
to a minimum of the energy in the kink sector, it is stable under perturbations which remain in that 
sector; in other words, under perturbations which preserve the kink boundary conditions on ip. The 
location of the kink is determined by the continuous parameter (zero-mode) b. The kink is highly 
localized, in the sense that its energy density (the summand in eqn 7) is concentrated on just a 
couple of links of the lattice. For example, ii b — 0, then the proportion of the energy contributed 
by the two links [—h, 0] and [0, h] is 4ft,/ (4 -I- ft^); if ft = 1 this proportion is 80%, and as ft —> 2 it 
tends to 100%. A diagram illustrating the kink solution ipix) for ft = 1 is presented in figure 1. 

It is instructive to look at the limit ft = 2. In this case, the solution of the Bogomol'nyi equation 
is supported at a single lattice site xq: 

if X < Xo 

ip{x) = c if X = xq (13) 
if X > Xo 

where c is a real constant. By contrast with (11), the "position" zero-mode is discrete; and there 
is an additional continuous zero-mode, namely the value c of V' at the centre of the kink. This 
limiting case is therefore marginally unstable: a perturbation will in general cause "0(3^0) to increase 
or decrease without bound. 
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Figure 1: The static kink for h — 1. 

The situation we wish to study is that of the moving kink. There is no expUcit solution in this 
case, and so one has to resort to approximation, or to numerical solution of the equations of motion 
(9), namely 

■ip = -^[sin(V'+ — ^) - sin(7/i — ^_)] - -[sin(7/'+ + ^) + sin(7/' -f ^_)] 
/i^ 4 
4 - /j2 4 + 

= cos V; (sin -I- sin •0-) sinV'(cosV'+ + cost/i-). (14) 

3 A collective-coordinate approximation 

The simplest approximation of kink dynamics is obtained by restricting tp to have the form (11), 
with b now becoming a dynamical variable b{t). So the number of degrees of freedom is reduced 
from infinity to one. Our conjecture is that this appproximation is a good one, provided that the 
speed V of the kink is small (cf [Q, ^); we shall elaborate on this condition later in the section. So 
the Lagrangian is 

L = Eji^ — Ep 

= f{b)b^-l, (15) 

where 

/(^) = V ^ sech\ix~b). (16) 

a; — — oo 

Note that f{b) is periodic with period /i, is even, and has a local maximum at 6 = 0. Graphs of f{b) 
for various values of h are given in figure 2. 

The Euler-Lagrange equation of the system is 

fib)b^-f'{b)b\ (17) 

This may be reduced to quadratures: 



= Fh{b), (18) 
where 6(0) = 0, 6(0) ~ v. The function Fh{b) is strictly increasing, and satisfies 

Fhib + nh) = Fhlb) + nFh{h) (19) 
for all integers n owing to the periodicity of /(6), so it sufhces to know Fh{h) for < 6 < /i. 
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Figure 2: The function f{b) of equation (1i 




Figure 3: The function b(t) over one wobble period for h = 1 and h = 1.8. 

This function F^ is easily inverted (fig. 3) to give the kink trajectory 

b{t)^F^\vt). (20) 

The main feature is that the kink wobbles as it moves through the lattice. This wobble is a dynamical 
effect; there is of course no static potential (such as Peierls-Nabarro) causing it. The time taken for 
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the kink to move from x — to x — h (the wobble period) is 

T=M^; (21) 

V 

and Fh{h)/h < 1, approaching the upper bound in the continuum hmit, /i — > 0. That is, the kink 
travels faster by a factor h/ Fh{h) than would be naively suggested by the initial velocity v. A graph 
of Fh{h) (fig. 4) shows that discreteness effects are small for h < 1.2, but grow large as h approaches 
2. 

We believe that this approximation of kink motion is accurate for small v. The Bogomol'nyi 
bound is crucial to this belief: the idea is that the configuration sits at the bottom of a one- 
dimensional potential valley, and the kink can move freely along the valley floor (ie. in the direction 
of the zero-mode b) . There are at least two ways in which the approximation differs from the full 
model. First, there is dynamical dressing (a moving kink has a different shape from a static one). 
In the continuum limit, this correction is of order v'^ (relativistic contraction), and one might expect 
the dressing in the discrete case also to be of order v"^ (bearing in mind that the expression for the 
static kink is exact). Secondly, there is transfer of energy to radiation (phonons). In Bogomol'nyi 
field theories, the suggestion is that this energy transfer will be suppressed by a factor of order 
exp(— Tq/Ti), where Tq is a typical timescale for the truncated dynamics, and Ti is the maximum 
period of the radiation in the theory Our suggestion (not proved) is that the same is true in 
Bogomol'nyi lattice models such as the present one. In this case, Tq is the wobble period (21), which 
is of order h/v; and Ti = 2tt, since the minimum phonon frequency is unity (see next section). So the 
energy transfer (per lattice site traverse) is suppressed by exp[— /i/(27rw)], and is therefore negligible 
for small v (as long as h is not small). If h is small, in fact if h is less than 1, then the amplitude of 
the wobble is tiny, and so one would expect the radiation to be negligible irrespective of v; of course, 
this is exactly what happens in the continuum (h 0) limit. 
1 f 




h 

Figure 4- Effect of discreteness on the function Fh{h)/h. 

The accuracy of the approximation has been tested numerically using a fully-explicit fourth-order 
Runge-Kutta algorithm with fixed time-step 0.01. The initial condition was a Galilean-boosted static 
kink profile with initial velocity 0.01 lattice sites per unit time [v — O.Ol/i). Simulations of duration 
1000 time units were performed for h ~ \, 1.2, 1.4, 1.6 and 1.8. In every case the kink moves freely, 
without pinning, undergoing motion of the predicted periodicity. Furthermore, inspection of the 
kink velocity over a single period reveals close agreement with h{t) calculated from (20) (see fig. 5), 
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although for h = 1.2, numerical errors rather swamp the very small theoretical wobble (note the 
scale on the velocity axis). 

Collective- variable analyses of the Frenkel-Kontorova model (cf have introduced the kink 
position as an extra variable, accompanied by a constraint. But since the dynamics is much more 
complicated, all the degrees of freedom have to be kept in order to obtain accurate results. A 
truncation to one degree of freedom, as in the present case, does not work. 
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Figure 5: Kink velocity over a single wobble period. 



4 Fast-moving kinks 

The collective-coordinate approximation is expected to fail at high velocities (except for small h). 
This is observed in the simulations as a gradual onset of kink deceleration as the initial velocity 
is increased. Failure occurs at lower velocities for coarser lattices — around v = 0.012 for h = 1.8 
compared with v — 0.15 for h — 1. The kink energy is dissipated in the form of small amplitude 
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oscillations ("radiation" or "phonons") emitted in its wake, propagating backwards (see fig. 6). 




Figure 6: Radiation by a fast-moving kink. 

The effect of this radiation on the kink velocity over a long time-scale can be seen in fig. 7. The 

data were produced by the above-mentioned Rungc-Kutta algorithm, run for 16000 time units with 
a time-step of 0.01. The initial configuration was a static kink Galilean-boosted to speed v = 0.3, on 
a lattice of unit spacing {h = 1). To cut reflexion of radiation from the fixed left-hand boundary, the 
first five lattice sites were damped. After an initial velocity drop of 0.02 in 10 time units as the kink 
assumes a more appropriate shape, it decelerates more slowly. The modulation of the amplitude of 
velocity oscillations is due to the velocity sampling (once every 10 time units) falling in and out of 
phase with the periodic wobble of the kink as it passes lattice sites. 
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Figure 7: Radiative kink deceleration. 
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The most interesting feature of fig. 7 is the existence of a threshold velocity, v « 0.16, below 
which deceleration, and hence radiation, is much reduced. Some understanding of this phenomenon 
may be gained by an analysis (motivated by Q) of the linearized equations of motion, 

•• 4 - ft,^ 4 + /i^ 

^--iX^(V^++^-)-^^- (22) 

From this one derives a dispersion relation for small-amplitude travelling waves, namely, 

2 4 + fe^ A-h^ 
" 2^^=°^^^' (23) 

where k is the wave-number. The angular frequency uj ranges between 1 and 2/h, the lower bound 
being responsible for the threshold velocity. 

As the kink travels along the lattice, it hits lattice sites with frequency v/h sites per unit time. 
Provided v > h/2'K, the kink can excite radiation of the same frequency, that is = 2ttv /h. However, 
if w < h/2'K, then w < 1, and the lattice cannot support such radiation. The kink can only excite 
higher harmonics, so the rate of energy dissipation is suddenly cut and the kink velocity becomes 
"quasi-stable." For h — \, this threshold occurs at w = 1/2tt « 0.159, in good agreement with the 
numerical data (fig. 7). 

A theoretical understanding of the specific shape of the graph is more elusive. We make the 
ad- hoc assumption (motivated by figure 7 of that the energy A_E lost by a kink in traversing a 
single lattice cell at speed v obeys an exponential law: 

= eP"-', (24) 

where p and q are positive constants, properties of the lattice. We further assume that the kinetic 
energy of the kink is 

Ek - ^cz;2 (25) 

where c is approximately constant provided h is not too large, as suggested by the results of section 
3 (fig. 2). These two equations imply a first order differential equation for v[t)^ easily solved to give 

v{t)^v{Q)- A\og{Bt+\), (26) 

where 

A ^ K 

P 

B = iLeP^W-?. 
he 

The dashed curve in fig. 7 is a fit of this formula to the numerical data, taking A = 0.032, B = 
0.004365, ^(O) — 0.28. Estimating c = 1.0986 by averaging the function f{h) for h — \, one deduces 
that p = 31.3, and q = 17.5. The fit is good for velocities greater than the radiation threshold at 
V = 0.159. 



5 Conclusion 

We have described a spatially-discrete sine-Gordon system, which is significantly different from the 
usual (Frenkel-Kontorova) system. There is a "topological" lower bound on the kink energy, and an 
explicit static kink solution with this energy. The kink moves on a "level playing-field" : there is no 
Peierls-Nabarro potential barrier. In particular, the behaviour of a kink moving with low speed is 
much simpler. 
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The idea described in section 2 for generating "topological" discrete systems, can be applied to 
many other continuum theories which have a Bogomol'nyi bound on the energy of a kink. By way 
of example, consider the 0"* equation, where the continuum potential energy density is 



1 ( dipV , 1 2\2 



In equation (4), one replaces cos^A by \^p^ — ^p- The most obvious factorization is then 



D = Aip, 



The lattice potential energy is then given by (5) as before, and the kink energy is bounded below 
by 2/3 (in the kink sector tp — > ±1 as — > ±cx)). 

It would be interesting to investigate breathers and kink-antikink collisions in this sine-Gordon 
model. Again, the absence of the Peierls potential should lead to behaviour which differs qualitatively 
from the usual model One simple observation is that in the limiting case h = 2, there exists a 
breather located at a single lattice site. Indeed, the field 
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